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KfiaHTOBHii MeTOfl BHSHaneHHa npoHHKHOCTi 

y FRW MO^ejii 3 pafliai^iero 

y po6omi poazAfidaem'bCH sanpuma Keanmoea 
ModeAb 0pidMaHHa-Po6epmcoHa-yoAKepa 3 nosu- 
muBHOK) KOCMOAOziuHOfo KOHcmaHmoH) , padiav^iew 
ma zasoM HariAuzma. /I^ah ananisa UMoeipHocmi 
HapodofceHHH ececeimy na iH^AHV^iuHiu cmadii hk 
(pyHKV,ii enepzii padiavfli yeedeno Hoee eusHaueuHH 
xeuAi, lUiO "eiAbHo" noiuupmembCH y cuAhHux no- 
ARX. SKoppenmoeaHO zpanuHuy yMoey myHeAwea- 
HHR BiACHKma, o6HUCAeHO npoxodMceHHR ma eid- 
6ummH y noemcmm KeanmoeoMy cmav,ioHapHOMy 
nidxodi. 

KAWHoei CAoea: (fiisuKa panubozo ececeimy, iH(pA}i- 
yiAepa-/le Bimma, zaa HanAuzina. 

E-mail: maidan@kinr.kiev.ua 

1 Introduction 

To date, among all variety of models of early Uni- 
verse one can select two prevailing approaches: 
these are the Feynman formalism of path integrals 
in multidimensional spacetime, developed by the 
Cambridge group and other researchers, called 
the "Hartle- Hawking method" [Ij, and a method 
based on direct consideration of tunneling in 4- 
dimensional Euclidian spacetime called the "Vi- 
lenkin method" [2j. In the quantum approach we 
have the following picture of the Universe creation: 
a closed Universe with a small size is formed from 
"nothing" (vacuum), where by the word "nothing" 
one refers to a quantum state without classical space 
and time. A wave function is used for a probabili- 
stic description of the creation of the Universe and 
such a process is connected with transition of a wave 
through an effective barrier. 

In majority of models tunneling is studied 
in details in the semiclassical approximations (for 
example, see [3H6]). Here, a tunneling boundary 
condition could seems to be natural, where the 
wave function should represent an outgoing wave at 
the large scale factor a. However, whether is such a 
wave free in asymptotic region? If to draw attenti- 
on on increasing gradient of potential, used with 
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opposite sign and having a sense of force, acting 
"through the barrier" on this wave, then one come 
to contradiction: influence of the potential on this 
wave is increased strongly at increasing a Now 
a new question has been appeared: what should the 
wave represent in the cosmological problem? So, we 
come to necessity to define "free" wave inside strong 
fields. 

The problem of correct definition of the wave 
in cosmology is reinforced more, if to calculate inci- 
dent and reflected waves before the barrier. Even wi- 
th known exact solution for the wave function there 
is uncertainty in determination of these waves. But 
penetrability is based on them. 



In order to estimate probability of formation of 
Universe accurately, we need in the fully quantum 
basis. Aims of this paper are: (1) to give definiti- 
on of the wave in strong fields; (2) to construct the 
fully quantum stationary method of determination 
of the penetrability and refiection using the defini- 
tion of the wave above; (3) to estimate how much 
the semiclassical approach is differed from the fully 
quantum one. 
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2 A model in the Friedmann Robertson 
Walker metric with radiation and 
generalized Chaplygin gas 

Let us start from of a closed {k = 1) FRW 

model in the presence of a positive cosmological 
constant A > 0, radiation and component of the 
Chaplygin gas. The minisuperspace Lagrangian has 
the form |8|: 



C (a, d) 



3 a 
8nG 



.0 , Svr G n 
-a + K a p[a^ 



(1) 



where a is scale factor, a is derivative of a with 
respect to time coordinate t, p (o) is a generalized 
energy density. In order to connect the stage of Uni- 
verse with dust matter and its another accelerating 
stage, in Ref. [9] a new scenario with the Chaplygin 
gas was applied to cosmology. A quantum FRW- 
model with the Chaplygin gas has been constructed 
on the basis of equation of state instead of p (a) = 
p{a)/3 (where p{a) is pressure) by pch = ~^/Pch' 
where A is positive constant and < a ^ 1. Soluti- 
on of equation of state gives 



PCh{a) 



A + 



B 



7,3 (1+a) 



1/(1+") 



(2) 



where B \s a, new constant. This model through one 
phase a connects the stage of Universe where dust 
dominates and DeSitter stage. At limit a — )• eq. ([5]) 
is transformed into pdust plus p\. From such limit 
we find A = p/i^, B = pdnst a-nd write the following 
generalized density |7]: 



p{a) 



PA + 



Pdust 
7,3 (l+«) 



1/(1+") 



+ 



Prad 



(3) 



where /5rad(fl) is component describing the radiation 
(equation of state for radiation is p{a) = /Orad(Q^)/3, 
p is pressure) and p\ = A/(87rG). 

The passage to the quantum description of 
the evolution of the Universe is obtained by the 
procedure of canonical quantization in the Dirac 
formalism for systems with constraints. We obtain 
the Wheeler-De Witt (WDW) equation, which after 
multiplication on factor and passage of component 
with radiation p^ad into right part transforms into 



the following form [7]: 

^ + ^ («)| V'(a) = ^rad ^{a), -Erad = 



Via) 



AttG 
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Pdust 
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(4) 
(5) 



a3(i+") 

where <^{a) is wave function of Universe. For the 
Universe of closed type (at k = 1) at SttG = = 
1 we have: 

Via) = 36a^-12a^{A + -^y''^"' (6) 

and E'rad = 12prad- 

Let us expand the potential ([6]) close to arbi- 
trary point a by powers q = a — a and restrict 
ourselves by linear item: 



Vch (q) = Vo + Viq. 
For coefficients Vq and Vi I find: 

Vo = Vbh(a = d), 
V = 72a + 12a3{-4A--^}: 



A + 



Pdust \-"/(l+") 



(7) 



(8) 



a3(i+") 
and eq. Q obtains form: 

- ^ ^(q) + (^0 - i^rad + Vi q) ^{q) = 0. (9) 
After change of variable 

Erad — Vo 



|yi|2/3 |yi|2/3 



we have 



(10) 



(11) 



3 Motivations to correct Vilenkin's 
boundary condition of tunneling 

Let us analyze how much a choice of the boundary 
condition in the asymptotic region is motivated. 

• In tasks of decay in nuclear and atomic physi- 
cs the potentials of interactions tend to zero 
in the asymptotic region. Here, an application 
of the boundary condition at limit of infini- 
ty does not give questions. In cosmology we 
deal with another, principally different type 
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of the potential: with increasing of the scale 
factor a modulus of this potential increases. 
A gradient of the potential, used with opposi- 
te sign and having a sense of force acting on 
the wave, increases also. So, there is nothing 
mutual with free propagation of the wave in 
the asymptotic region [7j. 

• Results of Ref. [SJ reinforce a seriousness of 
such a problem: the scale factor a in the 
external region is larger, the period of osci- 
llations of the wave function is smaller. This 
requires with increasing a to decrease step 
of calculations of the wave function. This 
increases time of calculations, increases errors. 
So, boundary condition in the asymptotic 
region has no practical sense in cosmology. 
In contrary, in nuclear and atomic physics 
calculations of the asymptotic wave are the 
most stable. 

• It has not been known whether Universe 
expands at extremely large scale factor a. 
Just the contrary, it would like to clarify this, 
imposing Universe to expand in initial stage. 

So, I redefine the boundary condition as |3| 

the boundary condition should fix the 
wave function so that it represents the 
wave, interaction between which and the 
potential barrier is minimal at such a 
value of the scale factor a where action 
of this potential is minimal. 

To give a mathematical formulation for this defini- 
tion, we are confronted with two questions: 

1) What should the free wave represent at arbi- 
trary point inside cosmological potential? 

2) At which coordinate is imposition of this 
boundary condition the most corrected? 

At first, let us solve the second question. Whi- 
ch should this point be: or this is a turning point 
(where the potential coincides with energy of radi- 
ation), or this is a coordinate where a gradient of 
the potential (having a sense of force of interacti- 
on) becomes zero, or this is a coordinate where the 
potential becomes zero? 



We define this coordinate where the force 
acting on the wave is minimal. We define 
the force as the gradient of the potential 
used with opposite sign. 

4 Definition of the wave minimally 
interacting v^ith the potential 

Definition 1 (strict definition of the wave). The 

wave is such a linear combination of two partial 
solutions of the wave function that the change of the 
modulus p of this wave function is closest to constant 
under variation of a: 



da^ 



p{a] 



(12) 



-dtp 



For some types of potentials it is more convenient 
to define the wave less strongly. 

Definition 2 (weak definition of wave). The wave 
is such a linear combination of two partial solutions 
of wave function that the modulus p changes mini- 
mally under variation of a: 



d_ 

da 



p{a] 



0. 



(13) 



a=atp 



Now we shall look for the function as 

^(e) = r.*(+)(C), (14) 

M/(±)(0 = J exp±z(-^ + /(On) dn, 



(15) 

where T is an unknown factor, /(^) is an unknown 
continuous function satisfying /(^) ^ const at — )> 
0, and Umax is the unknown upper limit of integrati- 
on. The real part of /(^) gives a contribution to the 
phase of the integrand function while the imaginary 
part of /(^) deforms modulus. 

At small enough values of |^| we represent /(^) 
in the form of a power series: 



(16) 



n=0 
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where /„ are constants. Substituting formula (fTC|) of the Numerov method with a constant step. By 



for ^'(^) into equation pT]) . we find (see Ref. [7J) 



p(±) 



p(±) 



Jn+2 



where 



2i 



X 



(±) 



(±)' 



4/1/2 '^2 ~ "Jq 



6i Ji 



(±) 



such a way, I obtain two partial solutions for the 
wave function and their derivatives in the whole 
studied region. 

Having obtained two linearly independent parti- 
al solutions 971(a) and ip2{a), we make up a general 
solution (prime is for derivative with respect to o): 



^ (n - m + l)(m + 1) /, 

m=0 



(±) f{±) 

n—m+l J m+1 



J. 



(±) 



i (n + l)(n + 2) 



(17) 



(±) 



J- 



(±) 



exp zbi 



u exp ibi 



(-y + /o'") c?-", (18) 



J, 



(±) 



(*ma: 

/ 



e 



du. 



In order to be the solution ^{i) closer to the well- 
known Airy functions, Ai (^) and Bi(,^), we choose 
/o = 0, /i = 1. 

5 Calculations of the v^ave function 

In order to provide a linear independence between 
two partial solutions for the wave function effecti- 
vely, I look for the first partial solution increasi- 
ng in the region of tunneling and the second one 
decreasing in this tunneling region. At first, I define 
each partial solution and its derivative at a selected 
starting point, and then I calculate them in the 
region close enough to this point using the method 
of beginning of the solution. Here, for the partial 
solution which increases in the barrier region, as the 
starting point I use arbitrary point a inside well with 
its possible shift at non-zero energy E^^^ or equals 
to zero a = at null energy -Brad, and for the second 
partial solution which decreases in the barrier regi- 
on, I select the starting point to be equal to external 
turning point atp, out- Then both partial solutions 
and their derivatives I calculate independently in 
the whole required range of a using the method of 
continuation of the solution, which is improvement 



if{a) = T- (Ci if I (a) + C2 if 2 (a)) , 



(±) ' Ci 



fHp'^ - f>\if2 



Co 



^Vi - 



(19) 



(20) 

where T is normalization factor, C\ and C2 are 
constants found from the boundary condition. 

6 Problem of interference betv^een the inci- 
dent and reflected waves 

Rewriting the wave function (/Jtotal in the internal 
region through a summation of incident (^inc wave 
and reflected c/^i-ef wave: 



V'total — V'inc + ^^ref ; 



(21) 



we consider the total flux: 

j(¥'total) = jinc + jref + Jmixed, (22) 

where 



(23) 



iinc = i Up\nc^ifX^c - h. C. j , 
jref = i (^V3refV(/J*gf - h. C. j , 

The 

Jmixed Component describes interference 
between the incident and reflected waves in the 
internal region (let us call it as mixed component 
of the total flux or simply flux of mixing). From 
constancy of the total flux jtotal we find flux jtr for 
the wave transmitted through the barrier, and: 

Jinc — Jtr Jref Jmixed, Jtr — Jtotal — COnst. 

(24) 

Now one can see that the mixed flux introduces 
ambiguity in determination of the penetrability and 
reflection for the same known wave function. 

In the radial problem of quantum decay definiti- 
on of penetrability and reflection looks to be condi- 
tional as the incident and reflected waves should be 
defined inside internal region from the left of the 
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Energy 


Penetrability Ppenctrability 


Time 


r 




Full QM method 1 


Full QM method 2 


Method WKB 


Full QM method 1 


Method WKB 


20.0 


7.6221 X 10^ 


-30 


7.7349 X lO^^y 


1.4692 X 10" 


-29 


9.4313 X lO+^y 


4.8928 X 10^ 


-29 


40.0 


3.5680 X 10- 


-26 


6.5169 X 10-25 


1.2298 X 10" 


-25 


2.1257 X 10+26 


6.1670 X 10^ 


-25 


60.0 


2.0591 X 10" 


-22 


4.3423 X 10-21 


8.1523 X 10- 


-22 


3.8719 X 10"*'22 


9.7797 X 10^ 


-21 


80.0 


1.5530 X 10- 


-18 


2.3850 X 10-^'^ 


4.4346 X 10- 


-18 


5.3862 X 10"^!^ 


1.8862 X 10^ 


-18 


100.0 


3.2922 X 10" 


-14 


1.1053 X 10-13 


2.0304 X 10- 


-14 


2.6622 X 10+1^ 


4.3167 X 10^ 


-14 


120.0 


8.6052 X 10- 


-11 


4.4005 X 10-1° 


7.9523 X 10" 


-11 


1.0678 X 10+11 


1.1555 X 10^ 


-11 


140.0 


2.2012 X 10 


-8 


1.5460 X 10-^ 


2.7128 X 10" 


-7 


4.3888 X 10+® 


3.5612 X 10" 


f7 


160.0 


2.9685 X 10" 


-5 


4.9980 X 10-3 


8.1663 X 10" 


-4 


3.4471 X 10+^ 


1.2530 X 10 


f4 


170.0 


3.4894 X 10" 


-3 


2.6078 X 10-1 


4.2919 X 10" 


-2 


3.0460 X 10+3 


2.4820 X 10" 


f2 



Ta6ji. 1: The penetrability Ppcnotrabiiity of the barrier and duration r of the formation of the Universe defined by 
eq. (|29|) in the FRW-model with the Chaplygin gas obtained in the fully quantum and semiclassical approaches 
(minimum of the hole is -93.579 and its coordinate is 1.6262, maximum of the barrier is 177.99 and its coordinate 
is 5.6866): the fully QM method 1 is calculations by the fully quantum approach for the boundary located in the 
coordinate of the minimum of the internal hole (i. e. coordinate is 1.6262), the fully QM method 2 is calculations 
by the fully quantum approach for the boundary located in the internal turning point atp. in (coordinates of the 
turning points are in Tabl. 2) 



barrier. In order to formulate these coefficients, we 
shall include into definitions coordinates where the 
fluxes are defined (denote them as xj^ft and Xj-ight)' 



right ) 



imc(a^left) 



R 



Jrcf(.a:icftj 

jmc(a;ieft) ' 



M 



J mixed (a^ieft) 
jinc(a;ieft) 



(25) 

From eqs. (|24|) and (|25p we obtain [7j {jtr and jj-ef 
are directed in opposite directions, jmc 

and jtr — in 

the same directions): 



ITI + \R\- M = 1. 



(26) 



Now we see that condition \T\ + \R\ = 1 has sense 
in quantum mechanics only if there are no any 
interference between incident and reflected waves 
which we calculate, and it is to use jmixed = 0. 

7 The penetrability and reflection: the fully 
quantum approach versus semiclassical 
one 

Now we shall estimate by the method described 
above the coefficients of penetrability and reflection 



for the potential barrier ([6]) with parameters A = 36, 
i? = 12 A, A = 0.01 at different values of the energy 
of radiation -E^ad- We shall compare the found coeffi- 
cient of penetrability with its value, which the semi- 
classical method gives. In the semiclassical approach 
we shall consider two following deflnitions of this 
coefficient: 



pWKB,(l) 
penetrability 



pWKB,(2) 
penetrability 



1 



(27) 



(29 + i/{2ey 



where 



(oxt) 



exp / \V{a) — E\ da. 



(28) 



(int) 



One can estimate also duration of a formation of the 
Universe, using by definition (15) in Ref. |10) : 



r = 2 atp, i 



1 



int 



penetrability 



(29) 
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Energy 


Fully quantum method 


Turning points 




Penetrability 


Reflection 


Interference 


Summation 


'^tp, in 


Otp, out 


20.0 


7.6221543404 


X 


10- 


-30 


1.00000000000000 


8.55 X 10" 


-2U 


1.00000000000000 


3.59 


7.05 


40.0 


3.5680158760 


X 


10' 


26 


1.00000000000000 


2.34 X 10" 


-19 


1.00000000000000 


3.79 


6.97 


60.0 


2.0591415452 


X 


10- 


-22 


1.00000000000000 


4.86 X 10- 


-20 


1.00000000000000 


3.98 


6.88 


80.0 


1.5530040238 


X 


10- 


-18 


1.00000000000000 


2.08 X 10- 


-19 


1.00000000000000 


4.18 


6.78 


100.0 


3.2922846164 


X 


10- 


-14 


0.99999999999996 


3.13 X 10" 


-20 


1.00000000000000 


4.38 


6.67 


120.0 


8.6052092530 


X 


10- 


-11 


0.99999999991394 


1.06 X 10" 


-19 


1.00000000000000 


4.59 


6.55 


140.0 


2.2012645564 


X 


10" 


-8 


0.99999997798735 


1.60 X 10- 


-19 


1.00000000000000 


4.83 


6.39 


160.0 


2.9685643504 


X 


10 


-5 


0.99997031435611 


6.94 X 10- 


-20 


0.99999999999961 


5.11 


6.18 


170.0 


3.4894544195 


X 


10 


-3 


0.99651054553176 


2.02 X 10" 


-19 


0.99999999995131 


5.31 


6.02 



Ta6ji. 2: The coefficients of the penetrability, reflection and mixing calculated by the fully quantum method and test 
on their summation for the FRW-model with the Chaplygin gas density component (the fully quantum approach 1 
is used at the internal boundary located in the coordinate of the minimum of the internal hole) 



Results are presented in Tabl. [TJ -P^r|^|^abiiity 
not included in the table because it coincidence wi- 
th -Pj^il^l^abiiity to the first 8 digits. One can see that 
the fully quantum approach gives the penetrability 
close to its semiclassical value, which differs from 
results [To] . 

In the next Tabl. 2 the coefficients of the 
penetrability, reflection and mixing calculated in the 



fully quantum method are presented for the energy 
of radiation -Erad close to the height of the barri- 
er. One can see that summation of all such values 
for coefficients allows to reconstruct the property 
(|26p with accuracy of the first 11-18 digits. Now 
it becomes clear that the approach proposed in 
Ref. |10| and the semiclassical methods do not gi- 
ve such an accuracy. 



8 Conclusions 

In the paper the closed Friedmann-Robertson- 
Walker model with quantization in the presence 
of a positive cosmological constant, radiation and 
Chaplygin gas is studied. Note the following. 

1) A fully quantum definition of the wave 
propagating inside strong field and interacti- 
ng with them minimally has been formulated. 
The tunneling boundary condition has been 
corrected. 

2) A quantum stationary method of determinati- 
on of penetrability and reflection relatively 
the barrier has been developed. Here, non-zero 
interference between the incident and reflected 
waves has been taken into account and for its 
estimation the coefficient of mixing has been 
introduced. 



In such an approach the penetrability of the barri- 
er for the studied FRW-model has been estimated. 
Note the following. 

• The probability for birth of asymptotically 
deSitter Universe is close to results obtai- 
ned by the semiclassical approach, but di- 
ffers on results obtained by non-stationary 
approach |10] (see Tabl. 1 and 2). 

• The refiection from the barrier has been 
determined at first time. It is differed essenti- 
ally on 1 at the energy of radiation close 
enough to the barrier height (see Tabl. 2). 

• The modulus of the coefficient of mixing (indi- 
cating interference between the incident and 
reflected waves) is less 10^^^. 

• A property (|26p is reconstructed to the first 
11-18 digits (see Tabl. 2). 
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